Inverse cascade in two-dimensional turbulence : deviations from Gaussianity 
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High resolution numerical simulations of stationary inverse energy cascade in two-dimensional turbu- 
lence are presented. Deviations from Gaussianity of velocity differences statistics are quantitatively 
investigated. The level of statistical convergence is pushed enough to permit reliable measurement of 
the asymmetries in the probability distribution functions of longitudinal increments and odd-order 
moments, which bring the signature of the inverse energy flux. Their scaling laws do not present 
any measurable intermittency correction. The seventh order skewness is found to increase by almost 
two orders of magnitude with respect to the third, thus becoming of order unity. 



PACS number(s) : 47.10. +g, 47.27.-i, 05.40. +j 

The inverse energy cascade in two dimensional Navier- 
Stokes turbulence is one of the most important phenom- 
ena in fluid dynamics. In agreement with the remark- 
able prediction of R.H. Kraichnan in 1967 M, the cou- 
pled constraints of energy and enstrophy conservation 
make the energy injected into the system flow toward 
the large scales. This is a basic difference with respect 
to 3D turbulence, where energy flows toward small scales 
in a direct cascade. The dynamical process of structur- 
ing and organization of the large scales by the inverse 
cascade is also of great interest for geophysical fluid dy- 
namics. First numerical and experimental observations 
of the inverse cascade and the ensuing Kolmogorov en- 
ergy spectrum were obtained in &K3|. The important 
point foreseen in [El is that the smallness of the skew- 
ness suggests that intermittency might be weak. This 
conjecture was later supported by numerical simulations 
p0| , pl| and experiments Q : scaling laws are compat- 
ible with dimensional predictions and both transversal 
and longitudinal velocity probability distribution func- 
tions (pdf's) look not far from Gaussian. The evidence 
stemming from experiments and simulations is that the 
inverse transfer takes place via clustering of small-scale 
equal sign vortices. Strong deviations from Gaussian- 
ity appear if the system has a finite size and friction 
extracting the energy from the large scales is small (or 
absent). A pile-up of energy akin to the Bose-Einstein 
condensation takes then place in the gravest mode ["J, 
large scale vortices are formed and energy spectra steeper 
than the Kolmogorov one are observed fl3|| . Here we 
shall not consider the condensation phase, concentrating 
on the inverse cascade statistics. Theoretically, inverse 
cascade enjoys a great advantage with respect to the di- 
rect one : the limit of molecular viscosity v — > can 
be taken without any harm in the equations of motion 
for velocity structure functions. At variance of 3D tur- 
bulence, the energy dissipation ^((Vv) ) is indeed van- 
ishing when v —* 0. The absence of dissipative anoma- 
lies is the clue for the analytical solution of inverse cas- 
cades in passive scalar advection 0^3], Intermittency 



was found to be absent, even though the statistics might 
be strongly differing from Gaussian. For 2D Navier- 
Stokes inverse cascade, dissipative terms can again be 
discarded but the situation is complicated by pressure 
gradients. They couple indeed the statistics of velocity 
differences S r v = v(r) — v(0) at various r's in a non- 
local way. Closures on velocity increments-pressure gra- 
dients correlations have been proposed by invoking the 
quasi-Gaussianity of the statistics and quantitative pre- 
dictions have been derived in this way jl(|. The issue 
of quasi-Gaussianity is however moot as deviations are 
intrinsically entangled to the dynamical process of in- 
verse energy cascade. Standard calculations (see, e.g., 
fll7|) permit indeed to derive the 3/2 Kolmogorov law 
for 2D turbulence: S L (r) = ([S r v ■ r] ) = 3/2er, where 
r = r/r. The energy flux is denoted by e and the fact 
that it goes upscale reflects into the positive sign of the 
moment. Precise quantitative informations on the devi- 
ations from Gaussianity are however difficult to obtain. 
Odd-order structure functions involve for example strong 
cancellations between negative and positive contributions 
and the 3/2 law itself could not be observed in previous 
studies, due to lack of resolution and/or statistical con- 
vergence. It is our purpose here to present the results 
of high-resolution numerical simulations aimed at quan- 
titatively analyzing deviations from Gaussianity in the 
inverse energy cascade. 

Specifically, the 2D Navier-Stokes equation for the vor- 
ticity oj(r, t) = —Aip(r, t) is : 



dtUJ + J (to, ip) = vIAlo — ato — A/, 



(1) 



where ip is the stream function, the velocity v = V^V' = 
(Vyip, —W x ip) and J denotes the Jacobian. The fric- 
tion linear term —ato extracts energy from the system at 
scales comparable to the friction scale rj{ r ~ e 1 / 2 ^ 3 / 2 , 
assuming a Kolmogorov scaling law for the velocity. To 
avoid Bose-Einstein condensation in the gravest mode 
we choose a to make rjf T sufficiently smaller than the box 
size. The other relevant length in the problem is the 
small-scale forcing correlation length If, bounding the 



inertial range for the inverse cascade as If <C r <C %■• 
We use a Gaussian forcing with correlation function 
(f(r,t)f(0,f)) = S(t-t')F(r/l f ). The ^-correlation 
in time ensures the exact control of the energy injection 
rate. The forcing space correlation should decay rapidly 
for r 3> If and we choose F(x) = F l 2 exp(—x 2 /2), 
where Fq is the energy input. The numerical integration 
of (pi) is performed by a standard 2/3-dealiased pseu- 
dospectral method on a doubly periodic square domain 
of iV 2 = 2048 2 grid points. The viscous term in (pi) has 
the role of removing enstrophy at scales smaller than If 
and, as customary, it is numerically more convenient to 
substitute it by a hyperviscous term (of order eight in our 
simulations). Time evolution is obtained by a standard 
second-order Adams-Bashforth scheme. After the sys- 
tem has reached stationarity, analysis is performed over 
twenty snapshots of the velocity field equally spaced by 
one large-eddy turnover time. 

Let us now discuss the results. In Fig. [Owe present the 
third-order longitudinal structure function S y L ' (r) com- 
pensated by the factor l/(er), showing a neat plateau at 
the value 3/2 - in agreement with the Kolmogorov law 
- over a range of almost one decade of scales. In Fig. g 
it is presented the energy spectrum E(k), which displays 
a clear Kolmogorov scaling k~ 5 > 3 , and the energy flux 

n(fc). 
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FIG. 1. Compensated third order longitudinal structure 
function S L (r)/(er). The dotted line is the value 3/2. Note 
the linear vertical scale. 



The Kolmogorov constant in 

E(k) = Ce 2/3 k- 5/3 , 



(2) 



is found to be C = 6.0 ± 0.4 . Previous numerical 
simulations and experiments report values of the Kol- 
mogorov constant C ranging from 5.8 to 7.0 MM P~3[ - 
The structure function constants corresponding to (0) 



(2) 



3C^ 2) /5 
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rC = 12.9 ± 0.8, where 



are C L - ju r 7 ^ — 2 s/3r(4/3) 2 

the first two equalities follow from isotropy and incom- 

pressibility and 
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FIG. 2. Energy spectrum E(k). The dashed line has a 
Kolmogorov power law scaling. In the inset is shown the 
energy flux II (k). 

For transverse moments, r is substituted in ([|) by r±, 
perpendicular to it. It is of interest to remark that lon- 
gitudinal and transverse velocity increments are uncor- 
related, i.e. S V L ' T (r) = 0. The relatively large value of 

(2") 

C L implies a small skewness of the longitudinal velocity 
differences (3/2)/(cf°) 3 / 2 = 0.03. Albeit the longitu- 
dinal pdf looks close to Gaussian and quite symmetric, 
nevertheless on a more quantitative ground asymmetries 
turn out to be quite strong as shown by the two curves 
of Sf\r) and S ( p{r) in Fig. |. 




FIG. 3. Longitudinal structure functions of order 5 (lower) 
and 7 (upper). Power laws with exponents 5/3 and 7/3 are 
drawn for comparison. 

First, we can observe that their scaling behavior is in 
agreement with Kolmogorov predictions, without signifi- 
cant anomalous corrections. Second, the constants are 



C 



(5) 



130 and C 



(7) 



skewness C 



(->) 



/(C (2) )5/2 



14000, giving for the hyper- 
- 0.22 and C { p /(C^) 7 ' 2 ~ 



4 n) (r) = <M-f] n > = 4 B V) 



n/3 



(3) 



1.8. The error bars can be estimated from r.m.s. fluc- 
tuations of compensated plots and for the seventh order 
(which is of course the most delicate) they amount to 
20%. The increase of the skewness by almost two orders 



of magnitude from the third to the seventh order is par- 
ticularly informative. Holder inequalities apply indeed 
to absolute moments and odd-order moments (without 
absolute values as in (ph) might a priori even become 
smaller when their order increases. Our main motivation 
was precisely to find out whether odd-order moments 
were decreasing or increasing with the order. The an- 
swer to this question shows that hyperskewness is defi- 
nitely not a "small parameter" to be used in perturbative 
schemes for the statistical properties of the inverse energy 
cascade. The predictions in pi , although based on a clo- 
sure explictly invoking small deviations from Gaussian- 
ity, turn out to be compatible with the numerical results. 
This indicates that the closure is likely to be more robust 
and "nonperturbative" than its derivation might suggest. 
Another striking evidence for the importance of the 
longitudinal pdf asymmetries is provided in Fig. |4| 




FIG. 4. Antisymmetric structure functions S_ (r) of order 
n = 4,6. The solid lines have slopes 4/3 and 2, respectively. 
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FIG. 5. Antisymmetric part of the longitudinal velocity in- 
crements pdf, rescaled according to Kolmogorov scaling, at 
separations r = 0.05 (□), r = 0.075 (x), r = 0.1 (+), lying 
in the inertial range of scales. In the inset, the same data in 
linear-logarithmic scale. 



We 



consider here the antisymmetric part of the 



pdf V (o~v L {r)) — V (— Sv L (r)) (shown in Fig. ||) an d 
calculate "antisymmetric structure functions" such as 
Si 4) = ^u A {T{u)~V{-u)) du. Both the fourth 
and the sixth moment have a clean Kolmogorov scal- 
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ing S K _'(r) = C w (erY L/0 , with C y _'/(C^) 

and Ci 6) /(C[ 2) ) 3 ~ 0.6. This indicates that the non- 
Gaussian antisymmetric part, although visually small, 
has imprinted all the relevant scaling informations on the 
inverse cascade. 

From the graph in the inset of Fig. g it can be appre- 
ciated that the asymmetric part of the pdf has relatively 
consistent tails, which are the cause for the large observed 
values of the hyperskewness of order 5 and 7. 

To investigate possible dependencies on the type of 
forcing we also performed (shorter) numerical simulations 
with an injection rate characterized by the spectral corre- 
lation function (f(k, t) f(k', t')) = 6(t - t')6(k + k')6(l - 
klf) [pi. At variance with the former choice, this forcing 
is limited to a narrow bandwidth in Fourier space but 
its spatial correlations decay rather slowly. Odd-order 
structure functions and the antisymmetric part of the 
pdf do not show any visible dependence on the details 
of the energy input. Conversely, the symmetric part of 
the pdf of velocity differences and even order moments 
are more sensitive. For the forcing limited to a shell of 
wavenumbers, both the (symmetrized) longitudinal and 
the transverse pdfs are visually indistinguishable from 
Gaussian, as shown in Fig. |6j. 
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FIG. 6. Left: symmetric part of the longitudinal velocity 
difference pdf. Right: pdf of transverse velocity differences. 
The forcing is restricted to a band of wavenumbers and the 
data refer to a single velocity field snapshot. Gaussian distri- 
butions are shown as solid lines. 

Deviations of kurtosis and hyperkurtosis from their 
Gaussian values are small and compatible with those pre- 
sented in p^| . For the forcing localized in physical space, 
the far tails of the pdf at scales 0{lf) tend to be broader. 
This tendency is due to the formation of small vortices of 
size comparable to If, which generate large velocity dif- 
ferences - especially transverse ones - across a distance 



of the order of their size. The effect becomes of course 
negligible at scales larger than If but it might affect the 
quality and the extension of the scaling region for even 



order structure functions |18|. 

In conclusion, we have presented quantitative evi- 
dences for deviations from Gaussianity of the velocity 
increment statistics in the inverse energy cascade. Odd- 
order structure functions display a clean power law scal- 
ing compatible with classical Kolmogorov predictions. 
Numerical prefactors in adimensionalized structure func- 
tions are expected to be universal with respect to the 
forcing statistics and have been measured up to the sev- 
enth order. Despite the small value of the skewness, de- 
viations have been shown to be quite strong and the hy- 
perskewness of seventh order to be of order unity. Asym- 
metries in longitudinal velocity statistics should therefore 
be incorporated and treated systematically in theoretical 
models for the inverse energy cascade. 
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